Linearly Polarized Gluons and the Higgs Transverse Momentum Distribution 
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We study how gluons carrying linear polarization inside an unpolarized hadron contribute to the 
transverse momentum distribution of Higgs bosons produced in hadronic collisions. They modify the 
distribution produced by unpolarized gluons in a characteristic way that could be used to determine 
whether the Higgs boson is a scalar or a pseudoscalar particle. 
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It is sometimes said that the LHC is a 'gluon collider', 
because at high energies the gluon density inside a proton 
becomes dominant over the quark densities. Higgs pro- 
duction, in particular, predominantly arises from gluon- 
gluon 'fusion' gg — > H through a triangular top quark 
loop. QCD corrections to this process have been cal- 
culated with increasing precision [1-7], making it well 
understood. It is not commonly known however that the 
LHC is actually also to some extent a polarized gluon col- 
lider, since gluons can in principle be linearly polarized 
inside an unpolarized proton. Their corresponding distri- 
bution, here denoted by hf 9 and first defined in Ref. [8] , 
requires the gluons to have a nonzero transverse momen- 
tum with respect to the parent hadron. It corresponds to 
an interference between +1 and —1 helicity gluon states 
that would be suppressed without transverse momentum. 

So far the function hf 9 has not been studied exper- 
imentally, and consequently nothing is known about its 
magnitude. Only a theoretical upper bound has been 
given [8, 9]. Recently, several ways of probing hf 9 have 
been put forward, namely in heavy quark pair or dijet 
production [9], or in photon pair production [10], where 
in all cases the transverse momentum of the pair is mea- 
sured. One way in which linearly polarized gluons can 
manifest themselves in these processes is through az- 
imuthal asymmetries. However, it was found that they 
can also generate a term in the cross section that is inde- 
pendent of azimuthal angle. This happens when two lin- 
early polarized gluons, one from each hadron, participate 
in the scattering. In this way they can also contribute to 
production of a scalar particle, such as a scalar or pseu- 
doscalar Higgs boson, when its transverse momentum q T 
is measured. It has in fact been shown [11, 12] that such a 
contribution is generated perturbatively. In other words, 
if at tree level gluons are taken to be unpolarized, at order 
a s they will become to some extent linearly polarized. In 
the transverse momentum distribution of spin-0 particles 
produced in proton-proton collisions this will give rise to 
an additional contribution at order a 2 s , because of the 



double helicity flip involved (see Fig. 1). While this may 
be expected to make only a relatively modest contribu- 
tion, the function 9 is of non-perturbative nature and 
is present at tree level already. Therefore, a significant 
influence of linearly polarized gluons on the distribution 
of the produced particle at low q T is not excluded. 




FIG. 1: Gluon helicities in the gg — > H squared amplitude for 
unpolarized (left) and linearly polarized production (right). 

In light of this, we will investigate in this letter how 
the distribution of linearly polarized gluons may affect 
the transverse momentum distribution of Higgs bosons 
for q T <C ran , where ran is the Higgs mass. We shall ob- 
serve that linearly polarized gluons may in fact provide 
a tool to uncover whether the Higgs boson is a scalar or 
a pseudoscalar particle. Thus far, relatively few sugges- 
tions to this end have been put forward for the LHC, typ- 
ically using azimuthal distributions, for example in Higgs 
+ jet pair production [13] or in r pair decays [14]. The 
suggestion we put forward here does not involve measure- 
ments of any angular distributions. Instead, we will show 
that linear polarization of gluons simply leads to a mod- 
ulation of the Higgs transverse momentum distribution 
that depends on the nature of the Higgs particle. 

Transverse momentum dependent distribution func- 
tions (TMDs) of gluons in an unpolarized hadron are 
defined through a matrix element of a correlator of the 
gluon field strengths F M!/ (0) and F ua (^), evaluated at 
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fixed light-front (LF) time £ + = £-n = 0, where n is 
a lightlike vector conjugate to the parent hadron's four- 
momentum P. Decomposing the gluon momentum as 
p = x P + p T + p~n, the correlator is given by [8] 



"p n a f d(£-P) d 2 j T ip< 
' (p-n) 2 J (2tt)3 6 

x(P\Tr[F^(0)F^(O]\P)\ LF 



- -— J n^ v f 9 - ( P ^ P t -4- n» v P ^ \ h ±9 \ C\ 



Withp^ = _ p 2 ; g gv = g^_ p » n v/ p . n _ n »pv/p. n ^ and 

M the proton mass. ff{x,p 2 ) represents the unpolar- 
ized gluon distribution and h 1 9 (x,p^) the distribution 
of linearly polarized gluons. In (1) we have omitted a 
Wilson line that renders the correlator gauge invariant. 
As any TMD, hf 9 will receive contributions from initial 
and/or final state interactions, which make the gauge 
link process-dependent. Therefore, despite the fact that 
it is T-even, h t 9 can receive non-universal contributions, 
and its extraction can be hampered for processes where 
factorization does not hold, such as dijet production in 
hadron-hadron collisions [15-17]. Higgs production, on 
the other hand, is expected to allow for TMD factoriza- 
tion, just like the Drcll-Yan process. A more detailed 
study of this remains to be carried out. 

The calculation of the Higgs production cross section 
in the TMD framework closely follows Refs. [15, 18]. The 
generic contribution by gg — > H reads 
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For now, we assume on-shcll production of the Higgs 
particle, with q and Eh its momentum and energy. 
P a and Pb are the momenta of the colliding protons, 



q 2 /(2P a(b) -q). To low- 



s' = (P a + Pb) 2 , and x a(b ) 
est order, the hard partonic amplitude M. is given by the 
well-known formula [4] for the gg — > H triangle diagram: 



M% = i2 1 / 4 G 1 J 2 a s m 2 H g^ A H (r) / (Sn) 
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for a scalar Standard Model (SM) Higgs boson H°, 
where we only consider top quarks in the triangle, and 
where Gp is the Fermi constant, a s the strong coupling 
constant, r = m H /(4m 2 ) with the top mass m t , and 
A h (t) = 2 (r + (r - 1) J(r)) /t 2 with 
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For a pseudoscalar Higgs boson A with a simple coupling 
g t (V^Gf) 1 / 2 ™- 2 75 to quarks [19] we have instead 



M% = i2 1 ^G 1 J 2 a s m 2 H 4" A a (t)/(8it), 
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where is the two-dimensional Levi-Civita tensor, and 
Aa(t) = gt2J(r)/T. As mentioned above, QCD correc- 
tions to these amplitudes have been calculated. However, 
since our goal is to study the effect of linearly polarized 
gluons whose distribution h^ 9 is anyway unknown, we 
limit ourselves to the lowest-order expressions (3) and 
(5). This leads to the following expressions for scalar 
and pseudoscalar Higgs production: 
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where the upper (lower) sign refers to the scalar (pseu- 
doscalar) case, and where we have the TMD convolution 

C[wff] = J d 2 p aT J d 2 p bT S 2 (p aT +p bT - q T ) 

^w(p aT1 p bT ) f(x a ,p 2 aT ) f(x b ,pl T ) , (7) 
with the transverse momentum weight 
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We emphasize the sign difference in the C[wh h\ 9 h^ 9 ] 
term in Eq. (6), which may offer an opportunity to de- 
termine the parity of the Higgs boson. The terms involv- 
ing hf 9 have the model-independent property (q^ a )hh = 
/ d 2 q T {q 2 T ) a C[w H h± 9 h± 9 } = 0fora = 0, 1. This fea- 
ture points towards a distinctive transverse momentum 
distribution of the h^ 9 h^ 9 term with a double node in 
q T . We note that since (l)hh = 0, linearly polarized 
gluons do not affect the q T -integrated cross section. 

In the following we estimate the possible size of the 
contribution by linearly polarized gluons to the Higgs 
production cross section at tree level. Although the func- 
tion hf 9 itself is unknown, a model-independent positiv- 
ity bound for it has been derived in Ref. [8]: 
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valid for all x and p T . The maximally possible effect will 
be generated when this bound is saturated. Models may 
also shed light on the size of hf 9 . In the simple pertur- 
bative quark target model of gluon TMDs of Ref. [20] the 
function h^ 9 is found to possess the same characteristic 
1/x increase as the distribution of unpolarized gluons /f , 
which suggests that linearly polarized gluons may be as 
relevant at small x as unpolarized ones. Another recent 
model calculation [21] shows saturation of the positivi- 
ty bound for h^ 9 in heavy nuclei in certain transverse 
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momentum regions (Weizsacker- Williams model) or even 
over the full momentum range (dipole model). This sug- 
gests that saturation of the positivity bound at least lo- 
cally in x or p 2 T might not be an unrealistic assumption. 

We follow a standard approach for TMDs in the litera- 
ture (see [22]) and assume a simple Gaussian dependence 
of the gluon TMDs on transverse momentum: 
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where G(x) is the collinear gluon distribution and the 
width (p 2 ) is assumed to be independent of x. The bound 
(9) is directly satisfied by the form 
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We choose r — 2/3. The left panel of Fig. 2 shows the 
p T -dependence of ff and hf 9 for two values of the Gaus- 
sian width: (p 2 T ) = lGeV 2 and (p 2 T ) = 7GeV 2 . The 
latter value may be more appropriate at Q = run, cf. the 
Gaussian fit to fi(x,p 2 ) evolved to Q = Mz of Ref. [23]. 
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FIG. 2: ie/t; Gaussian distributions for ff andh± 9 (divided 
by G(x)) as functions of p T for two different values of (p%) . 
Right: Resulting ratio R = C[wh hf 9 hf g ]/C[ff ff]. 

It is straightforward to compute the convolution inte- 
grals appearing in Eq. (6) analytically: 



C[/f/f] 



G(x a ) G(x b ) 



-exp 



C[wh h^ g h^ 9 



= G{x a )G{x b ) x 
367r(p 2 ) 
3(g 2 ) 2 



(Pi) 16(P 2 )< 



exp 



(12) 



(13) 



Their ratio R is shown in the right panel of Fig. 2. It is 
a measure of the relative size of the contribution by lin- 
early polarized gluons and shows the anticipated double 
node of C[w H h^ 9 h 1 9 ]. It is evident that at least within 
our simple model, linearly polarized gluons have a sizable 
effect on the Higgs q T distribution. We stress again that 



the effect enters scalar and pseudoscalar Higgs produc- 
tion with opposite sign. If the effect is at or near the level 
shown by Fig. 2, it should easily allow to determine the 
parity of the Higgs boson, provided a sufficiently fine scan 
in q T is possible in experiment, to resolve both nodes. 

So far we have considered only the production of an 
on-shell Higgs boson. In reality the Higgs will decay into 
some observed final state, and there will be background 
reactions contributing to this final state that are not re- 
lated to the Higgs. These backgrounds may themselves 
be sensitive to linearly polarized gluons. We will now 
briefly consider one example of this, the Higgs decay into 
a photon pair. We reserve a more detailed study of final 
states such as 7Z, ZZ or WW for a future publication. 

After production in gg — > H, the two-photon decay 
of a SM Higgs occurs through a top quark or IF-boson 
triangular loop. The decay of a pseudoscalar Higgs is 
instead described by physics beyond the SM and hence is 
model-dependent. There are often no tree-level couplings 
to IF-bosons in this case [24], so here we consider only 
the top quark coupling. For both a scalar or pseudoscalar 
Higgs, the lowest-order amplitude can be written as [25] 
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where a is the electromagnetic coupling and Th the Higgs 
decay width. For a scalar Higgs, = gT V &\ a \ b , whereas 
r^ v = \ a it^ S\ a \ b in the pseudoscalar case, with A a , Af, 
the photon helicities. Ah(t) and Aa{t) are given in 
Eqs. (3) and (5) with f = s/(4m 2 ), where s — (p a + 
Pi,) 2 ~ x a XbS for gluon momenta p a , p b . Finally, 



^h^ 77 (s) = W(t W ) + -N c A H {f) . 
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with t w = s/(4m2y) and W(r) = -(2t 2 + 3t + 3(2r - 
1)J(t))/t 2 describes the contribution by the W trian- 
gular loop. We assume J r J 4_ >77 (s) = ^^^^(f). In the 
following we consider a relatively light Higgs mass m# = 
120 GeV with a small total width T H ~ 5 x 10~ 3 GeV [26] . 

As is well known, an important QCD background to 
photon pair production at high energies is generated by 
99 ~^ 77 y i a a quark box [27]. This subprocess was 
studied recently in the context of TMD factorization in 
Ref. [10]. Using its results, we add the two lowest-order 
amplitudes describing the box diagram and the Higgs 
resonance and extract the azimuthal-angle independent 
cross section: 



da" 



F 1 C[f 9 f 9 ]+F 2 C[w H hi 9 hi 9 }. (16) 



Here q = q a + qb is the momentum of the photon pair. 
dVl = d(pd cos 8 denotes the solid angle element for each 
photon, with the angles <j>, 9 defined in the Collins-Soper 
frame [10]. Fi and F 2 are calculated functions of 9 and 
the pair mass Q = y/s that we will not give here. They 
depend on the box and Higgs amplitudes. 
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FIG. 3: Ratio F2/F1 as a function of pair mass squared in a 
region around mu = 120 GeV for various 



We find that the box contribution dominates the pro- 
cess except when the photon pair mass is close to the 
Higgs mass. Figure 3 shows the effect of the box-Higgs 
interference on the ratio F2 /Fi as a function of Q around 
niH = 120 GeV for a scalar or pseudoscalar Higgs. Away 
from Q = tub (by a few hundred MeV) we find F\ ^> F2 , 
such that the additional term from linearly polarized glu- 
ons contributes at most 10% to the cross section, but on 
average around 1% or less. However, near Q = mu where 
the Higgs contribution dominates, we find F\ rs ±i*2- 
The ratio of the second to first term in Eq. (16) then 
becomes approximately the ratio ±R of Fig. 2. Fig- 
ure 3 suggests that a distinction between a scalar and 
pseudoscalar Higgs is possible, if the experimental res- 
olution of the photon pair mass Q is sufficiently good. 
Higgs bosons in extensions of the SM which typically have 
larger widths, would required less fine Q-binning. Also, 
for heavier Higgs bosons other final states such as WW- 
or ZZ-production may allow for a better Q-resolution. 
In any case the Q-bin size around the Higgs mass is to 
be chosen as small as possible to maximize the effects 
caused by linearly polarized gluons. 

We conclude that the effect of linearly polarized glu- 
ons on the Higgs transverse momentum distribution can 
in principle be used to determine the parity of the Higgs 
boson, provided hf 9 is of sufficient size. Of course, it 
could turn out that hf 9 is in reality smaller than in our 
model or that it exhibits nodes in x or p T , complicating 
the analysis. Our results thus provide additional mo- 
tivation for experimental studies of h^ 9 using different 
probes, such as dijet and heavy quark or photon pair 
production. We stress that perturbative gluon-radiation 
effects will alter the q T distributions expected on the ba- 
sis of our simple Gaussian model. Their inclusion will 
require merging our model with the soft-gluon resumma- 
tion techniques described in [11, 12, 23, 28, 29]. This 
will also affect the eventual size of the contribution by 
linearly polarized gluons. A full study of this is needed. 
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